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ABSTRACT 



Context. Diffuse radio emission, in the form of radio halos and relics, traces regions in clusters with shocks or turbulence, probably 
produced by cluster mergers. The shocks and turbulence are important for the total energetics and detailed temperature distribution 
within the intracluster medium (ICM). Only a small fraction of clusters exhibit diffuse radio emission, whereas a large majority of 
well-studied clusters shows clear substructure in the ICM. Some models of diffuse radio emission in clusters indicate that virtu- 
ally all clusters should contain diffuse radio sources with a steep spectrum. External accretion shocks associated with filamentary 
structures of galaxies could also accelerate electrons to relativistic energies and hence produce diffuse synchrotron emitting regions. 
The detection of radio emission from such filaments is important for our understanding of the origin of the Warm-Hot Intergalactic 
Medium (WHIM), and relativistic electrons and magnetic fields in the cosmic web. Here we report on Giant Metrewave Radio 
Telescope (GMRT) observations of a sample of steep spectrum sources from the 74 MHz VLSS survey. These sources are diffuse on 
scales > 15'', and not clearly associated with nearby (z < 0.1) galaxies. 

Aims. The main aim of the observations is to search for diffuse radio emission associated with galaxy clusters or the cosmic web. 
Methods. We have carried out GMRT 610 MHz continuum observations of unidentified diffuse steep spectrum sources. 
Results. We have constructed a sample of diffuse steep spectrum sources, selected from the 74 MHz VLSS survey. We identified eight 
diffuse radio sources probably all located in clusters. We found five radio relics, one cluster with a giant radio halo and a radio relic, 
and one radio mini-halo. The giant radio halo has the highest radio power (Pi. 4) known to date. By complementing our observations 
with measurements from the literature we find correlations between the physical size of relics and the spectral index, in the sense that 
smaller relics have steeper spectra. Furthermore, larger relics are mostly located in the outskirts of clusters while smaller relics are 
located closer to the cluster center. 
Conclusions. 
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1. Introduction 

Studies of large-scale structure (LSS) formation have made sig- 
nificant advances during the last decade. It has been found that 
nearly all massive clusters have undergone at least several merg- 
ers in their history and that presently clusters are still in the 
process of accreting matter. A significant fraction of the accret- 
ing mass is in the form of (smaller) clusters and galaxy groups. 
Cluster mergers are the most energetic events in the present day 
Universe, with kinetic energies of the order of 10^^ - 10^^ erg, 
which are dissipated in giant shock waves and turbulence. An 
important aspect is the total energy budget and the detailed tem- 
perature distribution within the ICM, both of them are a ffected 
by the merger history of a cluster (e.g. JPolag et al.ll2008h . 

Diffuse steep spectrum radio emission is observed in about 
50 massiv e merging and pos t-merging galaxy clusters (see the 
review by iFerrari et alJ [20081 and references therein). This dif- 
fuse emission is difficult to detect due to its low surface bright- 
ness and steep spectral inde43, so the fraction of clusters hosting 
diffuse radio emission is probably larger than we currently know. 
The diffuse emission i n clusters is com monly divided into three 
main classes (Feretti & Giovannini 1996). Radio Halos are ex- 
tended (> 1 Mpc) diffuse unpolarized (< 5%) sources, located 
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in the center of clusters. They have a regular smooth appear- 
ance, and follow the thermal X-ray emission. Radio Relics are 
elongated structures with an irregular morphology, mostly lo- 
cated in the periphery of clusters. Relics can be highly polarized 
(10 - 50%). S everal different subclasses have been identified 
( Kem pner et al. 2004). Most known radio relics and halos are 
found in clusters which show signs of a current or recent merger. 
This supports the scenario in which the relativistic electrons are 
accelerated by merger-induced shocks or turbulence. However, 
Radio Mini-halos are not associated with merging clust ers. They 
are found in the centers of coo l core clusters (e.g., Fa bian et aH 
I1991I: iPeterson & Fabianll2006h and are associated with the cen- 
tral cluster galaxy and typically have sizes < 500 kpc, with 
the diffuse emissio n surrounding the central cluster galaxy (e.g., 
iGovoni et ani2009h . 

Two different mechanisms for in- situ acceleration of parti- 
cles have been proposed to explain relics in clusters: (i) adiabatic 
compression of fossil radio plasma by a passin g shock wave 
producing a so called radio ''pho enix" (EnBlin & Gopal-Krishnal 
I2OOII: lEnBlin & Briiggenl l2002h . or (ii) diffusi ve shock ac- 
celeration (DSA) by the Fermi-I process (e.g., Drurv "1983'; 
Blandford & Eichler 1987; Jones & Elfison 1991; Ensslin et al] 
19981; iMalkov & O'C Drurvll200ll) . In the first scenario, radio 
relics should have toroidal and complex filamentary morpholo- 
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gies. These relics are capable of producing very steep, curved ra- 
dio spectra due to inverse Compton (IC) and synchrotron losses. 
In the DSA scenario the electrons are accelerated by multiple 
crossings of the shock front (in a first order Fermi process). 
These relics have large sizes (Mpc) and are direct tracers of 
shock fronts in clusters. The spectral index is determined by the 
balance between the continuous acceleration at the shock front 
and energy losses in the post-shock region. 

The diff'use emission within clusters reveals the presence of 
relativistic electrons and magnetic fields on scales ~ 1 Mpc. 
Spectral aging, due to synchrotron and IC losses of the emit- 
ting electrons may explain the steep spectra. If the electrons 
are injected via Fermi acceleration (DSA), their energy fol- 
lows a power-law distribution. The power-law index of the in- 
jecte d electrons is related to the Mach number of the shocks 
fe.g. JHoeft & Bruggenl l2007h: shocks with a low Mach number 
have steeper radio spectra. Clearly, low-freque ncy surveys are 
needed to l ocate and stud y these sources (.Cassano et al.ll2QQ6[ 
l2QQ7ll2Q08h . Interestinglv. lBrunetti et al.l (l20Q8h discovered a ra- 
dio halo in the cluster Abell 521, which was previously known 

to host a radio relic, with a spectral index of 2.1, suggesting 

the existence of a population of diff'use source in clusters with 
spectral indices < -1.5. 

Numerical simulations show the dev elopment of variou s 
types of shocks during structure formation (Mini ati et al.l l2QQQ). 
These shocks diff'er in their location with resp ect to t he clus- 
ter center and Mach numbers 4M iniati e t aP I2QQQI: iMiniatil 
l2QQ2t iRvu et al.ll2QQ3l: iPfrommer et al..2006t ' lVazza et al.l l2009')T 
External accretion shocks have A1 » 1 and process the low- 
density, unshocked intergalactic medium (IGM). This results in 
relatively flat spectral indices of about -0.5 at the location of the 
shock front. Further away from the shock front the spectral in- 
dex steepens due to synchrotron and IC losses. Internal shocks, 
(i.e., merger and flow shocks) occur within the cluster. The Mach 
numbers of these shocks are lower resulting in steeper spectral 
indices. Binary merger shocks are the result of a cluster merging 
with a another cluster or a larg e sub-structure. In this case a dou- 
ble radio relic is expected (e.g. jRoettiger et al.ll999l:lHoeft et al.l 
12008: van Weeren et al. 2009^^ 

As pointed out by ICen & Ostrikerl (Il999l) . hydrodynamic 
models indicate that up to half of the baryons at present 
time should have temperatures in the range of 10^ - 10^ K. 
Unfortunately, studying the abundance and distribution of this 
WHIM is very challenging, since its main tracers are highly 
excited Oxygen lin es which are difficult to observe (e.g., 
iNicastro et al.ll2005b . A fraction of the accretion shocks will be 
supersonic and can accel erate energetic electrons up to energies 
of 10^^ - 10^9 eV (e.g.. Nor man et al.lll995l: Ikang et al.lll996l: 
llnoue et al ]|200i). In the presence of magnetic fields, such elec- 
trons will emit faint diff'use synchrotron radiation. The detection 
of these radio filaments is very important as this would provide 
a probe of the WHIM. Recent magnetohydrodynamical mod- 
eling indicates that detecting radio emissio n from the filamen- 
tary cosmic web should be possible (e.g., iKeshet et"al l2004l: 
Hoeft & Brii ggen 2007|). iHoeft et al.l (l2008h : iPfrommeii (l2008l) : 



Pfrommer et al.i (2008) however find that in the outskirts of clus- 



ters (at a few times the virial radius) or filaments, external ac- 
cretion shocks cause little radio emission, owing to the low den- 
sity o f both magnetic field energy and cosmic ray (CR) particles 
there (iMiniati et al.ll200ll) . They are therefore difl&cult to detect 
even with the sensitivity of upcoming radio telescopes such as 
LOFAR. Relic emission from internal accretion shocks occur 
in a higher density environment so that they should be detected 
with current radio facilities. 



When searching for radio halos, relics and filaments in low- 
frequency radio surveys, various other steep spectrum sources 
are also present. These include ultra-steep (angular size < 15'') 
spectrum sources (USS, see IMiley & De Breuck 2008, for a re- 
view) associat ed with high-z radio ga laxies (HzRG), "fossil" or 
"dying" FR-I dFanaroff^cfe Rilevlll974 radio sources, and "head- 
tail" sources, the last two having a steep spectrum due to spectral 
aging of the radio emission. In high-resolution (< 5'') observa- 
tions, for example at 1.4 GHz with the Very Large Array (VLA), 
diff'use objects will be resolved out due to missing short base- 
lines. This provides a method for selecting diff'use radio sources 
associated with galaxy clusters or the cosmic web. FR-I sources 
can be partly excluded by removing sources that are clearly as- 
sociated with individual galaxies. 

The 74 MHz VLA low-frequency Sky Survey (VLSS), 
ICohen et aP (l2007l) , covers about 3n steradians of sky north of 
S = -30°. The resolution of the survey is 80'' (FWHM) and the 
rms noise level is about 0.1 Jy beam"^. The source catalog con- 
tains roughly 70, 000 sources with a point source d etection limit 
of 0.7 Jy beam-i. A new calibration algorithm (ICotton et al.l 
12004 was used to remove the ionospheric distortions, which 
can be severe at this low-frequency. The 1 .4 GHz NRAO VLA 
Sky Survey (NVSS^. IComion et al.l (Il998h . covers the entire sky 
above S = -40°. The NVSS images have a rms noise of about 
0.45 mJy beam"\ and a resolution of 45''(FWHM). The catalog 
contains about 2x10^ sources above a flux of ~ 2.5 mJy. 

In this paper we present radio continuum observations of 
26 difl'use (angular size > 15'') steep spectrum sources selected 
from the VLSS survey with the GMRT at 610 MHz. The main 
aim of this project is to determine the morphology of the sources 
and search for difl'use structures which could be associated with 
shock fronts or turbulence in clusters, and accretion shocks onto 
filaments of galaxies. 

The layout of this paper is as follows. In Sect. [2] we discuss 
the sample selection, this is followed by an overview of the ob- 
servations and data reduction in Sect. O In Sect. |4] we present 
the radio maps of the most interesting sources and discuss these 
sources individually. By combining our radio observations with 
data from the literature (X-ray and optical observations) we have 
tried to classify the sources. In Sect. [51 spectral indices are mod- 
eled using our flux measurements combined with literature val- 
ues. We end with a discussion and conclusions in Sects. [6] and [71 

Throughout this paper we assume a ACDM cosmology with 
ifo = 71 km s-i Mpc-\ = 0.3, and Qa = 0.7. 

2. Sample Selection 

Spectral indices (between 1400 and 74 MHz) were calculated 
for all sources in the VLSS survey. We selected sources with 
a < -1.35 which resulted in a total of 176 sources. This cutoff" 
is somewhat arbitrary, but a significant lower cutoff" resulted in 
a very small number of sources selected while a higher cutoff 
would select too many sources for follow-up observations. VLA 
B-array 1.4 GHz ~ 5 min snapshot observations were carried 
out on March 25-29 and May 10, 2005 of a subsample of 68 
from the 176 sources. Two intermediate frequencies (IFs) with 
a bandwidth of 50 MHz each were used, centered at 1385 and 
1465 MHz. From this 68 sources, 36 were found to be resolved 
out. This showed that these sources had extended emission on 
scales > 15''. From the 36 sources 13 were identified with 
known nearby galaxies. The remaining 23 sources were included 
in the sample. We have also searched for additional sources 
with a < -1 .15, by making use of the 1.4 GHz FIRST survey 
(5'' FWHM, iBecker et al.lll995h . Sources with a FIRST ffux at 
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least 8 times lower than the 1.4 GHz NVSS flux were initially 
selected and visually inspected to remove obvious double lobe 
sources. The spectral index cutofl' of -1.15 was chosen because 
higher values resulted in too many double lobes to be selected 
for visual inspection. Furthermore, most known radio halos and 
relics have spectral indices steeper than this value. The amount 
of flux resolved out in the selection criterium was a tradeofl' as 
lower values also resulted in too many sources to be selected for 
visual inspection. After visual inspection we found three addi- 
tional sources which showed the presence of difl'use emission. 
None of these sources were clearly associated wit h nearby indi- 
vidu al galaxies in the POSS-II or SDSS surveys ( Abazajian et al. 
1200 9). The final list of sources and their coordinates are given in 
Tableffl 

3. Observations & Data Reduction 

High- sensitivity radio observations at 610 MHz were carried out 
with the GMRT in February and November 2008 of a sample 
of 26 diff'use steep spectrum radio sources. We divided a total 
102 hours of observation time evenly between the 26 sources. 
A total of 32 MHz bandwidth was recorded, using both the up- 
per (USB) and lower sidebands (LSB) which included both RR 
and LL polarizations. The data were collected in spectral line 
mode with 128 channels per sideband (IF), resulting in a spec- 
tral resolution of 125 kHz per channel. To increase UV-coverage 
we cycled between various sources, typically spending 40 min 
on a sources before moving to the next source. However, due to 
scheduling constraints it was only possible to do this for about 
half of our sources. The observations resulted in a net on-source 
time of ~ 3 hours, after flagging certain time-ranges which were 
afl'ected by radio frequency interference (RFI) or had other prob- 
lems. 

The data were reduced and analyzed with the NRAO 
Astronomical Image Processing System (AIPS) package. 
Bandpass calibration was carried out using the standard flux cal- 
ibrators: 3C48, 3C147, and 3C286. Fluxes of 29.43 Jy (3C48), 
38.26 Jy (3C147), and 21.07 Jy (3C28 6) at 610 MHz were as- 
signed to th ese sources using t he Perley & Taylo r (1999) exten- 
sion to the iBaars et aP (Il977h scale. A set of 6 channels free 
of RFI was taken to normalize the bandpass (channel 15 - 20) 
for each antenna. Strong RFI was removed automatically (with 
the AIPS task 'FLGIT'). The data was then visually inspected 
for remaining low-level RFI using the AIPS tasks 'SPFLG' and 
'TVFLG' . After that an initial phase and amplitude calibration 
was carried out using the bandpass and secondary calibrators, 
where we also transferred the flux densities from the primary cal- 
ibrators to the secondary calibrators. The found solutions were 
then transferred to the target sources. We have not chosen to av- 
erage any channels in order to minimize the efl'ects of bandwidth 
smearing and to aid further removal of RFI. The first and last few 
channels of the data were discarded as they were noisy. ^ 

Fo r making imag es v^e use d the polyhedron method (Perley 
119891: ICornwell & Perleylll992t) to minimize the eff'ects of non- 
coplanar baselines. Both USB and LSB were simultaneously 
gridded, imaged and cleaned. We used a total of 199 facets to 
cover ~ 2 times the full primary beam. This made the removal 
of sidelobes from strong sources far away from the field center 
possible. After a first round of imaging, in some cases "ripples" 
were seen in the maps which were subsequently removed af- 
ter identifying the corresponding baseline(s). Several rounds of 
phase self-calibration were carried out before doing a final am- 
plitude and phase selfcalibrati on. Images were made using ro- 
bust weighting (robust = 0.5. iBriggslll995l) and corrected for the 



primary beam attenuation. Images were cleaned to 3 times the 
rms noise level to minimize clean bias eff'ects. 

The thermal noise in each map is expectecfl to be 

CTthermal - ^ , , ,,,, , = , UJ 

G ^|n{n - l)NiFAytint 

with = 92 K the system temperature, G = 0.32 K Jy"^ the 
antenna gain, ^ 28 the number of working antennas, Nif = 2 
the number of sidebands used (both recording RR and LL po- 
larizations). Ay = 13.5 MHz the bandwidth per sideband, and 
^int the net integration time. The expected thermal noise for 3 hrs 
integration time is about 32 yuJy beam"\ where we have taken 
into account that typically 20% of the data is flagged due to RFI. 
The noise levels in our maps range from 40 to 202 yuJy beam"^. 
The noise levels dependent on the UV-coverage and the pres- 
ence of strong confusing sources in the field of view limiting 
the dynamic range. The lowest noise level of 40 yuJy beam"^ 
is quite close to the thermal noise level. The uncertainty in the 
calibratio n of the ab solute flux-scale is between 5 - 10%, see 
Chandra eTaP (120041) . 

4. Results 

All sources, except VLSS J2229. 1-0136, were detected in the 
GMRT images. The NVSS counterpart of VLSS J2229.1-0136 
has a position ofl'set of about 40'', the flux density in the NVSS 
and VLSS surveys of these sources is also close to the detec- 
tion threshold. We therefore conclude that VLSS J2229. 1-0136 
is most likely a noise peak as the GMRT observations should 
have easily detected the source. From here on we have left out 
VLSS J2229.1-0136 in the further discussions. A summary of 
the beam parameters and noise levels for the maps is given in 
Table [T] To identify optical counterparts overlays were made us- 
ing SDSS and POSS-II images. Spectroscopic redshifts were in- 
cluded from the literature. For sources without a spectroscopic 
redshift, but having SDSS DR7 coverage we took the SDSS pho- 
tometric redshift. For other sources we used that K and R mag- 
nitudes of massive elliptical galaxies correlate with redshift (K -z 
and R-z relations, e.g., Will ott et al.ll2003l:lde Vries et al.ll2007h . 

We have also checked for any X-ray counterparts to the radio 
sources. In the next subsection we describe the most interesting 
sources that appear to be related to shocks or turbulence in clus- 
ters or filaments. The other sources are discussed in AppendixlAl 

4.1. Individual sources 
4.1.1. VLSS J1 133.7+2324 

The radio emission shows two parallel filamentary structures and 
a northern patch of diff'use emission connected with the west- 
ern filamentary component. Some faint radio emission is extend- 
ing towards the west. This extension coincides with the galaxy 
UGC 6544 (MCG+04-27-065; PGC 35694) located a z=0.02385 
(Havnes et al. 1997|). This galaxy was clas sified as a Sbc spiral 
by Nilson (1973) and iPaturel et al.l (l2003l) . The source is also 
Hsted in the IRAF Faint S ource Catalog (Fl 131 1+2341) and 
detected at 60 and 100 yum (iMoshir & et al.lll990l) . The galaxy 
has been included in 21 cm hydrogen emission line surveys. 
Springob et al. (2005) give a self-absorption corrected line flux 
of 2.97 + 0.32 Jy km s"^ and a line- width of 372 + 3 km s"^ .The 

^ http://www.gnirt.ncra.tifr.res.in/gmrt_hpage/ 
Users/doc/manual/UsersManual/node 1 3 .html 
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Table 1. Source properties and results 



Source Name 


RA 

(J2000) 


DEC 
(J2000) 


Beam size 
arc sec 


^rms 

^Jy 


TT 

^610 

mJy 


^ 1400MHz 
"^74MHz 


LAS^^ 
arcsec 


Redshift 
z 


VLSS J0004.9-3457 


00 04 53.63 


-34 56 34.1 


11.4x6.6 


73 


158 + 16 


-1.40 + 0.04 


116 


0.29 + 0.08^ 


VLSS J0227.4-1642 


02 27 26.74 


-16 42 47.2 


5.3 X 4.0 


94 


37.0 + 6.2 


-1.50 + 0.06 


89 


>0.7 


VLSS J0250.5-1247 


02 50 30.83 


-12 47 30.3 


9.4 X 3.7 


118 


33.2 + 3.3 


-1.94 + 0.07 


21 


>0.7 


VLSS JOS 11.6+0254 


05 11 37.79 


+02 54 19.4 


8.7x6.1 


139 


98.0+10.3 


-1.42 + 0.05 


45 


0.20 + 0.05^ 


VLSS J0516.2+0103 


05 16 17.89 


+01 03 40.1 


8.1 x6.5 


79 


17.4 + 4.4 


-1.73 + 0.06 


35 


>0.7 


VLSS J0646.8-0722 


06 46 52.10 


-07 22 37.9 


5.7 X 4.8 


63 


120+13 


-1.70 + 0.05 


59 


0.23 + 0.06^ 


VLSS J0717.5+3745'' 


07 17 30.92 


+37 45 29.7 


8.2 X 6.0 


78 


501 + 50 


-1.15 + 0.04 


171 


0.5548^ 


VLSS J0915.7+25ir 


09 15 41.51 


+25 11 48.2 


8.6x5.9 


134 


194 + 21 


-1.52 + 0.04 


63 


0.324" 


VLSS Jl 117. 1+7003 


11 17 06.46 


+70 03 57.1 


7.8x4.3 


52 


9.1+0.9 


-1.87 + 0.07 


18 


>0.7 


VLSS Jl 133.7+2324 


11 33 44.73 


+23 24 50.6 


6.6 X 3.9 


51 


80.6 + 10.1 


-1.69 + 0.06 


84 


0.61+0.16^ 


VLSS J1431.8+1331 


14 31 49.90 


+ 13 31 54.1 


5.3 X 4.8 


51 


120+13 


-2.03 + 0.05 


101 


0.159936" 


VLSS J15 15. 1+0424' 


15 15 09.18 


+04 24 41.1 


7.6 X 5.4 


79 


192 + 33 


-1.50 + 0.05 


169 


0.0972^ 


VLSS J1636.5+3326 


16 36 34.93 


+33 26 33.9 


5.2x4.7 


40 


16.1 + 1.9 


-1.84 + 0.14 


28 


0.65 + 0.19^ 


VLSS J1710.5+6844 


17 10 35.03 


+68 44 58.4 


8.9 X 4.7 


63 


23.3 + 2.8 


-1.75 + 0.06 


45 


0.28 + 0.08^ 


VLSS J1930.4+1048 


19 30 27.19 


+ 10 48 02.5 


5.3 X 4.8 


56 


67.4 + 7.0 


-1.90 + 0.06 


101 




VLSS J2043.9-1118 


20 43 58.46 


-11 18 45.6 


5.8 X 4.2 


77 


53.9 + 6.3 


-1.74 + 0.05 


41 


0.43 + 0.15^ 


VLSS J2044.7+0447 


20 44 43.64 


+04 47 24.5 


10.5 X 7.5 


202 


57.7 + 5.8 


-1.55 + 0.06 


71 


0.46 + 0.15^ 


VLSS J2122.9+0012 


21 22 54.14 


+00 12 03.4 


5.8x3.8 


81 


24.0 + 2.4 


-2.00 + 0.08 


12 


> 1.4 


VLSS J2209.5+1546 


22 09 32.91 


+ 15 46 29.9 


6.9 X 6.2 


64 


36.4 + 3.8 


-1.56 + 0.07 


61 


>0.7 


VLSS J2213.2+3411 


22 13 12.45 


+34 11 51.6 


5.5 X 4.6 


94 


127 + 13 


-1.55 + 0.04 


62 


1.6 + 0.5^'^ 


VLSS J2217.5+5943 


22 17 30.39 


+59 43 05.3 


6.3 X 4.3 


47 


79.9 + 8.0 


-2.20 + 0.06 


99 




VLSS J2229.1-0136 


22 29 11.95 


-01 36 58.8 


8.6x4.2 


101 










VLSS J2241.3-1626 


22 41 22.57 


-16 25 35.7 


6.1 x5.8 


129 


69.7 + 7.3 


-1.44 + 0.06 


47 


>0.7 


VLSS J2341. 1 + 1231 


23 41 06.91 


+ 12 31 36.9 


5.5 X 4.4 


50 


117 + 14 


-1.70 + 0.04 


115 


0.62 + 0.15^ 


VLSS J2345.2+2157 


23 45 15.47 


+21 57 55.1 


4.8 X 4.2 


63 


79.3 + 8.4 


-2.14 + 0.05 


70 


0.23 + 0.06^ 


VLSS J2357.0+0441 


23 57 05.54 


+04 41 14.7 


7.8 X 4.9 


82 











" redshift estimated using the fitted Hubble-K relation from Willott et al. ( 2003) 
^ redshift estimated using the fitted Hubble-R relation from .de Vries et al.. (.20071) 
" spectroscopic redshift from SDSS DR7 

^ median photometric redshift from SDSS DR7 

^ using the band magnitude fromlPe Breuck et aDdloolB) 

^lEdge et al.' (^2003 |) 

^[struble&Rood ( 199^ 

^ LAS = largest angular size 

' resolved out in the 1.4 GHz FIRST survey 

^ fluxes are extracted from the regions indicated in the figures by dotted lines. 

^ for the compact sources (LAS < 45'') the fluxes were measured by fitting single Gaussians to the sources after convolving the maps to the 
45''NVSS resolution 



galaxy also has a small UV-excess (iMiyauchi-Isobe & Maeharal 
l2QQ2h . 

Using the IRAS fluxes we can estimat e the star forming rate 
(SFR) in this galaxy ( Solo mon et al.l 1997b by calculating the Far 
Infrared Luminosity (Lfir) 

^FIR = 3.94 X 10^(2.585' 60//m + S lOO//m)K^ 60fim/S 100//m)^L ' (^) 

with Lfir in Lq, the flux S in Jy, the luminosity dis- 
tance in Mpc, and ^5 60/^100) a color correction factor 
(iLonsdale&Heloul 119851) which is of the order of one. This 
gives a SFR of ~ 3Moyr- \ using SFR (Moyr'^) = Lfir/(5.8 x 
lO^Lo), iKennicuttI (Il998l) . The radio flux of the galaxy is dif- 
ficult to measure as part of the radio emission from the galaxy 
overlaps with the emission from the filamentary source. For the 
part which does not overlap with the steep spectrum source we 
measure a flux of about 6 mJy. We adopt a total flux for the 
galaxy of 10 mJy, assuming we missed about 40% of the pre- 
viously reported flux. The luminosity at 1.49 GHz (Li.49ghz) is 
L = 47rDlSy(l + zy~\ with Sy the observed flux at the rest 
frequency, and a 0.5. This results in L1.49GHZ = 8.2 x 10^^ 



W Hz"^ Using the FIR radio correlation (ICondon et al.lll99lh 
we find a luminosity of ~ 10^^ W Hz"\ in good agreement. 

The galaxy is unlikely to be associated with the filamen- 
tary radio structure to the east, given the above calculation, the 
steep spectral index, morphology, and spatial oflTset. In the back- 
ground (partly behind the spiral galaxy) there is an overdensity 
of faint red galaxies. These have a median photometric redshift 
of about 0.6 (SDSS DR7) and follow roughly the filamentary 
radio source. This is probably a cluster or a galaxy filament, 
with the foreground galaxy hiding part of the cluster/filament. 
Based on the morphology, the radio source is then classified as a 
relic. If the steep spectrum radio source is located at a distance 
of z = 0.6, the LAS of 84'' corresponds to a physical size of of 
660 kpc. 

4.1.2. VLSS J1 431 .8+1 331 

The radio emission shows two distinct components, the brightest 
one located to the east. This source has an optical counterpart, 
the cD galaxy of the cluster MaxBCG J217.95869+13.53470 
(iKoester et al.l l2007l) . The cluster is also detected in X- 
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Fig. 1. GMRT 610 MHz radio map. The Gaussian beam size is 
shown at the bottom left corner. Contour levels are drawn at 
V[l,2,4,6,8, 16, 32,...] x Acr^r^^, see Tabled The dotted lines 
indicate the region where the 610 MHz flux (85 lo), reported in 
Table [H has been extracted. 
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Fig. 3. GMRT 610 MHz radio map. Contour levels are drawn as 
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Fig. 2. SDSS r band image overlaid wit h GMRT 610 MHz co n- 
tours. Contour levels are drawn at V[l,8,32, 128,512,...] x 
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rays by ROSAT (IVoges et al.l 1 19991 l2000h as source 
IRXS J143150.6+133256. A SDSS DR7 spectrum puts 
the cD galaxy at a distance of z = 0.160. Three other sources 
in the cluster have a measured redshifts of 0.15962, 0.159153, 
and 0.164794. At this redshift the radio emission corresponds 
to a physical size of 170 kpc for the eastern and 125 kpc for 
the western component. The eastern component is an irregular 
curved structure with several bright knots. One of these cor- 
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Fig. 4. SDSS r band image overlaid with GMRT 610 MHz con- 
tours. Contour levels are drawn as in Fig.[2l 



responds with the nucleus of the cD galaxy. The morphology 
suggests that we are either seeing the interaction of radio plasma 
from the central AGN with the surrounding ICM or a small 
central radio source with a bright relic tracing a shock front 
(seen in projection). 

The western fainter component does not seem to be asso- 
ciated with any galaxy. The source could either be remnant ra- 
dio emission from a previous AGN episode or the signature of 
a shock wave. In the former case, the spectral index should be 
steeper then eastern component because of spectral aging. 
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Fig. 5. GMRT 610 MHz radio map. Contour levels are drawn as 
in Fig.O 



4.1 .3. VLSS J221 7.5+5943, 24P73 

This source was found to have an ultra- steep spectrum (a = 
-2.58 ± 0.14) by the Synthesis Telescope of the Dominion 
Radio Observat ory (DRAQ) Galactic plane survey at 408 MHz 
and 1.42 GHz (iHiggsl 119891: iJoncas & Higgslll990h. Since this 
source was located in the galactic plane Jjo ncas & Hig gs* ('1990') 
suggested the sources might be a pulsar. However, they noted 
that the source seemed to be slightly extended at 408 MHz. 
At 1.42 GHz the source was resolved into two distinct com- 
ponents. Both of these components were also listed as a single 
source 25P23. The sourc e was also d etected in the 38 MHz 8C 
survey (lReesl[T990: Hale s et al.lll995h . Subsequent L-band (1.4 
GHz) and X-band (8.4 GHz) observations with the VLA by 
iGreen & JoncasI (Il994l) . separated the source clearly into two 
components. A compact northern component and a southern 
diffuse component. The X-band observations only detected the 
northern component. By comparing the fluxes of the VLA and 
DRAG observations they concluded that the southern difl'use 
component had an ultra- steep spectrum and provided the bulk of 
the emission at low frequencies. Since the source was resolved 
this ruled out a pulsar identification. In fact it was suggested that 
this source might be a radio halo or relic located in a galaxy 
cluster behind the galactic plane. 

GMRT observations of 24P73 show a complex filamentary 
source. The source has so me similarities with the relic sources in 
AbeU 85 and Abefl 133 (Slee et al.ll2001 ). POSS-II and 2MASS 
images covering the area do not show the presence of any cluster. 
However this is not u nexpected given the extinction of = 6.7 
(ISchlegel et al.lll998l) . Given the steep spectral index and mor- 
phology we conclude that the sources is a radio relic. Deep NIR 
imaging will be necessary to identify the galaxy cluster associ- 
ated with the relic. 



4.1.4. VLSS J0004.9-3457 

The central component of the radio source is associated with 
a K= 14.86 mag elliptical galaxy. An overdensity of galax- 
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Fig. 6. GMRT 610 MHz radio map. Contour levels are drawn as 
in Fig.[T] 
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Fig. 7. POSS-II red image overlaid with GMRT 610 MHz con- 
tours. Contour levels are drawn as in Fig.O 



ies indicates the pre sence of a cluster, candidate B02291 
(IZanichelli et al.|[2QQl]) from radio-optically selected clusters of 
galaxies. We estimate a redshift for the cluster of 0.29 ± 0.08 
(using the K-band magnitude and the K-z relation). The cluster 
may be linked to a larger cluster about 3' towards the southwest 
at a redshift of 0.33 ± 0.09. The X-ray source 2E 0002.2-3515 
(iHarris et al.l 11996) is probably associated with this cluster of 
galaxies visible in POSS-II images. The diff'use radio emission 
surrounding the central elliptical galaxy resembles a mini-halo. 
No X-ray emission from the ICM of the cluster is detected in 
the ROSAT All Sky Survey. This implies that the cluster is not 
very massive. The object is similar to the mini-halo in the clus- 
ter MRC 0116+111 located at z = 0.131 (iGopal-Krishna et al.1 
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Fig. 8. GMRT 610 MHz radio map. Contour levels are drawn as 
in Fig.[Tl 



l2QQ2t iBagchi et ani2QQ9l) . which was also not detected in the 
ROSAT All Sky Survey. The size of th e mini-halo is abou t 
200 kpc, similar to other mini-halos (e.g. jGovoni et al ][200i. 
The radio image shows a bright knot at RA 00^ 04"^ 52.5' 
DEC -34° 56' 55'', to the south of the main component. This 
knot has an optical counterpart. To the east an arc-like struc- 
ture extends from the central component and bends towards the 
south. The arc is not associated with an optical counterpart. The 
origin of the arc is unclear, it could be part of a disturbed FR-I 
source, or a relic-like structure of fossil radio plasma. 

4.1 .5. VLSS J071 7.5+3745, MACS J071 7.5+3745 

The radio source is associated with the massive X-ray luminous 
cluster MACS J0717.5+3 745 at z^O.5548, with an overall ICM 
temperature of 11.6 keV (lEbeling et al.ll200ll: lEdge et al.ll2003l: 
lEbeling e t al .1120071). The cluster shows a p ronounced substruc- 
ture in optical images. lEbeling et aP (l2004l) reported the discov- 
ery of a large-scale filamentary structure of galaxies connected 
to the cluster. NVSS, WENSS and VLSS images reveal the pres- 
ence of a steep- spectrum radio source (a = -1.15). The ra- 
dio source was classified as a radio relic by "Edge et al.l (^2003). 
iMa et al.l (l2009l) presented X-ray (Chandra) and optical obser- 
vations (Hubble Space Telescope, ACS; Keck-II, DEIMOS) of 
the cluster. They found the cluster to be an active triple merger. 
Temperatures in the cluster exceeding 20 keV were found in 
some regions. Regions with a lower temperature of ~ 5 keV 
were found at the position of two subclusters, probably remnants 
of cool cores. 

Here we shortly describe our results for this cluster, a 
more detailed analysis in combination with additional archival 
VLA observations has been pre sented in a separate paper 
(Ivan Weeren et al.ll2009bh . see also Bonafede et al.l (l2009a ). Our 
radio maps reveal a complex source, consisting of diff'erent com- 
ponents. The main component is a twisted structure, with en- 
hanced regions of radio emission in the north and southwest 
of the cluster. These regions are connected by two bridges of 



Fig. 9. Chandra X-ray map overlaid with radio contours at 
610 MHz from the GMRT. The color scale represents the X-ray 
emission from 0.5 - 7.0 keV. The image has been adaptively 
smoothed using the TARA*^ package using a minimal signifi- 
cance of 5. Contour levels are drawn at [1, 2, 4, 8, 16, 32, ...] x 

0.312 mJy beam"^.^ 

http://www.astro.psu.edu/xray/docs/TARAI 
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Fig. 10. HST F814W ACS image overlaid with GMRT 610 MHz 
contours. Contour levels are drawn as in Fig.O 



emission to a bright central elongated source. The structure has 
a linear size of 700 kpc. No obvious counterparts are visible 
for both the north and south-west components. The presence of 
the two radio bridges suggest a relation with the central compo- 
nent. Although, no obvious counterpart was found for the cen- 
tral component by lEdge et al.l (|2003|) . we identify an elliptical 
galaxy at RA 07^ 17"^ 35^5, DEC -f37° 45' 0575 as a possi- 
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Fig. 11. GMRT 610 MHz radio map. Contour levels are drawn Fig. 12. SDSS r band image overlaid with GMRT 610 MHz con- 
as in Fig. [U tours. Contour levels are drawn as in Fig. O 



ble counterpart. If this is indeed the case then the source could 
be a wide-angle tail (WAT) source with the two bridges being 
the tails of the central source and the north and south-west com- 
ponents the hotspots. The eastern boundary of the structure is 
sharp, while on the other side some faint emission is seen ex- 
tending further westwards. Diffuse radio emission is also seen to 
the south of the main structure. The diffuse emission within the 
cluster has a total size of about 1 .2 Mpc and given that it roughly 
follows the X-ray emission we classify it as a radio halo. Clearly, 
the emission is not associated with individual sources. Using a 
spectral index of -1.2, typical for radio halos, we estimate the 
radio power (Pia) to be 5 x 10^^ W Hz"^. This makes it the 
most powerful radio known to date, in agreement with the X-ray 
luminosity-radio power (Lx - Pia) and temperature-radio power 
(T - P\ a) correlations (iLiang et al.ll2000l: lEnBlin & RottgeringI 
l2002HCassano et 2il\\20m . 

Towards the south a fainter linear structure is seen. A com- 
pact core, located halfway the linear structure is associated 
with a bright elliptical foreground galaxy (RA 07^ 17"^ 37^2, 
DEC -h37° 44' 23''). The source is probably a FR-I source asso- 
ciated with the compact core. 

Interestingly, the main twisted radio structure is located in 
between the brightest X-ray emission of the cluster. The main 
cool core visible in the X-ray image has no radio emission asso- 
ciated with it. The central radio structure also coincides with re- 
gions having a significantly higher X-ray temperature > 15 keV. 
The ICM temperature and X-ray morphology of the cluster fa- 
vor of a relic-like interpretation. We therefore conclude that the 
twisted radio structure is a giant relic tracing a shock front linked 
to the merger activity of the system. 

4.1.6. VLSS J091 5.7+2511 

The radio map shows a diffuse region of emission associated 
with a cluster of galaxies (MaxBCG J138.91895+25. 19876) at a 
redshift of 0.324. The radio source consists of a northern compo- 
nent and a fainter southern one. To the west a source is associated 



with an elliptical galaxy (RA 09^ 15^ 39^7, DEC +25° IV 37''). 
A few possible counterparts are visible in SDSS DR7 images for 
the northern diffuse component. The southern diffuse component 
has no obvious optical counterpart. This could be a radio relic 
with a projected size of about 190 kpc. High-resolution observa- 
tions will be needed to confirm this interpretation. 

4.1 .7. VLSS J1 51 5.1+0424, Abell 2048 

The radio image shows a filamentary radio source in the periph- 
ery of the cluster Abell 2048 at a redshift of 0.0972. The source 
consists of three elongated structures orientated roughly east- 
west. The three structures connect to the east. The source could 
be a complex double WAT source. However, the radio structure 
itself does not seem to be connected to any particular galaxy. The 
source has a projected size of 310 kpc, if located in the cluster. 
In the southeast a compact double-lobe source is associated with 
a large elliptical galaxy (RA 15^ 15^ 14M, DEC +04° 23' 10") 
located in the cluster. On the other side of the cluster a 0.19 Jy 
source (PMN Jl 5 15+0421) limits the dynamic rang e. The clus- 
ter is also detected in X-rays as RX J15 15.2+0421 (' Bade et al.1 
fT998h . A substructure on the east-side of the main cluster is vis- 
ible, hinting at a possible cluster merger, see Fig. [151 Given the 
location of the radio source at the edge of the cluster, the lack 
of a connection with a single optical counterpart, the steep ra- 
dio spectrum, and the morphology, we classify the source as a 
peripheral radio relic. 



5. Spectral index modeling 

Radio spectra can be an important tool to understand the origin 
of the relativistic electrons and to determine the age of the ra- 
dio emitting plasma. We have combined our flux measurements 
at 610 MHz with literature values to determine the radio spec- 
tra for the sources in our sample. We have included flux mea- 
surements from t he following surveys: 38 MHz 8C (Rees 1990; 
iHales et al.|[T995b . 74 MHz VLSS, 151 MHz 7C (.Waldram et all 
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Fig. 13. GMRT 610 MHz radio map. Contour levels are drawn 
as in Fig.[T] 
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Fig. 14. SDSS r band image overlaid with GMRT 610 MHz con- 
tours. Contour levels are drawn as in Fig.[2l 



Il996t iHales et al.ll20()7h. 232 MHz Mivun dzhang et al.lll997h . 
325 MHz WENS S ('RengelinketalJ 11997), 35 2 MHz WISH 
(iDe Breuck et al.1 l2Q02a) . 365 MHz TEXAS (iDouglas et al.l 
ll996|) and 1400MHz NVSS. 

We model the integrated radio s p ectra using the Jaffe- 
Perola (JP) model (Jaffe & Perolal Il973h descri bed by 
iKomissarov & Gubanovl (Il994|) . see also ISlee et iD (12001). 
We assume that (i) synchrotron self- absorption is negligible as 



Fig. 15. X-ray emission from ROS AT in the 0. 1 - 2.0 keV energy 
band. The image has been convolved with a circular Gaussian 
of 225''. The solid contours represent the radio emission at 
610 MHz from the GMRT. The radio contours are drawn at 
[1, 2, 4, 8, 16, 32, ...] X 7(Trms- Dashed contours show the galaxy 
distribution from SDSS DR7. Only galaxies with a photomet- 
ric redshift between 0.06 + Zerr < z < 0.15 - Zerr were selected 
from the catalogs, with Zerr the error in the photometric redshift 
and z = 0.0972 The galaxy isodensity contours are drawn at 
[6, 9, 12, 15, ...] galaxies arcmin"^. 



this only occurs in compact sources, (ii) radiative energy losses 
dominate over other losses such as adiabatic ones, (iii) the 
magnetic field is spatially uniform and constant in time, (iv) ra- 
diative electrons do not escape from the region, (v) the emission 
in the region is uniform, (iv) relativistic electrons were injected 
at a single point with a power-law distribution of energy, and 
(vii) the pitch angles of the synchrotron emitting electrons are 
assumed to be continuously isotropized on a timescale shorter 
than the radiative timescale. Relativistic electrons lose their 
energy by synchrotron emission and IC scattering off' the cosmic 

microwave background (CMB). 

In the Kardashev-Pacholczyk (KP) model ( Kardashev 1962'; 
iPacholczyk 1 970), also used by Ko missarov & Gubanov (1994), 
the pitch angle of the electrons remains in its original orientation 
with respect to the magnetic field. This introduces one more free 
parameter in the spectral index modeling (the ratio between the 
source magnetic field (B) and the eff'ective magnetic field for IC 
losses (^ic). We have chosen not to fit this model to limit the 
number of free parameters. Furthermore, the JP model is more 
realistic from a physical point of view, as an anisotropic pitch an- 
gle distribution will become more isotropic due to changes in the 
magnetic field strength between di fferent regions and scattering 
by sel f-induced Alfven waves (e.g. JCarim et al.ll 19911: ISlee et al.l 
1200 Ih . 

Our adopted scenario is as follows: When the source starts, 
we assume that it is fueled at a constant rate for a certain time 
^ci (the time of the continuous injection (CI) of relativistic elec- 
trons, with an injection spectral index c^inj). This is followed by a 
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Table 2. JP model fits with ain\ 
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= -0.75 



Source Name 




^CI 


to 








10^ yr 


10^ yr 


//Gauss 


VLSS J0646.8-0722 




0.227 + 0.0775 


0.0869 + 0.00624 


24.5 


VLSSJl 133.7+2324 




0.121 +0.034 


0.0870 + 0.00634 


23.3 


VLSS J1710.5+6844 




0.421 + 0.259 


0.0946 + 0.00812 


23.1 


24P73 (VLSS J2217.5+5943), z 


= 0.05 


0.156 + 0.0464 


0.120 + 0.00406 


29.1 


24P73 (VLSS J2217.5+5943), z 


= 0.10 


0.173 + 0.0534 


0.131 +0.00503 


27.0 


24P73 (VLSS J2217.5+5943), z 


= 0.15 


0.176 + 0.0557 


0.131 +0.00443 


26.4 


24P73 (VLSS J2217.5+5943), z 


= 0.20 


0.169 + 0.0500 


0.129 + 0.00497 


26.3 



relic phase (RE) where the injection of electrons is switched off' 
Ore)- During both of these phases the electrons lose energy by 
synchrotron and IC losses. The total source age is ^age = to + ^ci- 
In the spectrum two break frequencies occur both related to 
spectral aging. The first is the break frequency Vb of the first 
injected electron population (at the beginning of the CI phase) 



B 



(3) 



A second higher frequency break y/b is caused by the last elec- 
tron population injected, at the end of the CI phase 



y/b = Vb 1 + — 



(4) 



The injection spectral index ( ofinj) is determined by the AGN 
in the case of (relic) radio lobes or radio phoenices. For relics 
where the particles are accelerated by DSA the injection spectral 
index is related to the Mach number (At) of the shock. In linear 
theory the relevant expression is 



see Ivan Weeren et aP ('2009aV The injection time must be com- 
parable with the crossing time of the relic region by the shock 
front. Since "DSA-relics" are mostly located in the cluster pe- 
riphery, where the pressure of the ICM gas is lower, adiabatic 
energy losses may become important. These adiabatic energy 
losses aff'ect the radio spectrum, depending on the expansion rate 
of the relic and evolution of th e magnetic field strength (B(t)) 
with time (Kardashej Il96i iGoldshmidt & Rephaelil 1 19941: 
iMurgia et al. 2002). However, we only have a limited number 
of ffux measurements for our sources and only one proposed 
DSA-relic with enough ffux measurements to model the spec- 
trum. Furthermore, the location of this relic in the cluster as well 
as the identiff cation of the cluster itself are uncertain. We there- 
fore chose to ignore the eff'ects adiabatic expansion losses. 

The JP model is thus characterized by four free parameters: 
(1) the injection spectral index Ofinj , (2) the length of the CI phase 
tci, (3) the length of the RE phase to, and (4) a ffux normaliza- 
tion constant. To reduce the number of free parameters we hav e 
chosen to keep ofinj fixed to a value of -0.75 ('Parma e t"aDl2007h . 

The spectra are fitted by minimizing the squared value of 
the fit in a two-step process. We first determine the shape of the 
spectrum for diff'erent values of to and tci, both ranging from 
to 10^ yrs in 25 equal logarithmically spaced steps. Then an over- 
all ffux scaling (normalization) constant is determined by mul- 
tiplying the spectrum with a constant until the squared value 
is minimized. In this way, a 2-dimensional (25 x 25) array of 
squared values is created. We continue the fitting using the same 



process but now centering the to and tci values around the mini- 
mum ;^-squared value in the array and increasing the time resolu- 
tion by a factor of 1.5. The process is repeated until both to and 
^ci converge to a constant value, i.e., do not change by more than 
1% between subsequent iterations. The formal errors in the fit- 
ted parameters are determined by the corresponding distribution 
of the squared values. For the sources where we successfully 
fitted the JP-model we should take in mind that some of the sim- 
plifying assumptions we made may not be valid and could have 
aff'ected our results. 

For the magnetic field strength we use the revised 
equipartition magnetic field strength B' (iBrunetti et al.lll997l: 



Beck & Krause 2005). We use the same procedure as in 
van Weeren et al. ( 2009a) and take for the depth of the source 
(d) the average of the major and minor axis. The ratio of energy 
in relativistic protons to that in electrons (k) is set to 100. For the 
low-energy cutoff' (ymin, the energy boundary indicated by the 
Lorentz factor) we take 100. For other values of k and ymin, ^eq 

scales with y^^^''^^^^~''\l + y^)^/^^"^\ with a the spectral index. 

We have only fftted the radio spectra for sources with a red- 
shift (because the energy loss rate due to IC scattering is propor- 
tional to (I -\- z)^) and at least four ffux measurements available. 
An exception is 24P73 (or VLSS J2217.5+5943) for which we 
have no re dshift. If 24P73 h as a similar size as the relic in AS 5 
(150 kpc, ISlee et al.l[2QQlh . then its redshift would be around 
0.1. Since this is a very rough estimate we have also fitted the 
radio spectrum using redshifts of 0.05, 0.15, and 0.2. 

The ffuxes for the sources could be contaminated by the pres- 
ence of ffeld sources within the beam of the low-resolution sur- 
veys. In the case of 24P73, there is a 5.0 mJy NVSS source 
(NVSS J221736+594403) nearby which is blended with the dif- 
fuse source in the 8C, VLSS, and WENSS surveys. We fit- 
ted a second order polynomial in log-log space to the 610 
MHz GMRT, NVSS, and the 1.4 and 8.4 GHz ffuxes from 
iGreen & JoncasI 11994) for this source. We then extrapolated the 
ffux to 325, 74, and 38 MHz and subtracted off' this ffux to re- 
cover the uncontaminated ffux for the diff'use source. For the 
three other diff'use sources no unrelated sources were found that 
could have significantly contributed to the ffux. 

The fitted radio spectra are show in Figs. [161 to [191 The val- 
ues for tch to, and B^^ are reported in Table O The fitting pro- 
cess for the sources VLSS J1636.5+3326, VLSS J2213.2+3411, 
VLSS J0004.9-3457 did not converge. For VLSS J1636.5+3326 
the 1.4 GHz ffux measurement is relatively high, in this case 
the source may have restarted its activity causing the high 
1.4 GHz ffux value. This is also seen for several radio sources 
by Parma et al. ( 2007). For the other sources the ffux measure- 
ments are too closely spaced in frequency to provide enough 
constraints for the fitting process. 
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Fig. 16. Jaffe-Perola fit to the flux measurements. The duration Fig. 18. Same as Fig.[T6]but for VLSS J1710.5+6844. 
of the CI and RE phases are indicated in the figure. 
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Fig. 17. Same as Fig.[l6]but for VLSS Jl 133.7+2324. 
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Fig. 19. Same as Fig.IHbut for 24P73 (VLSS J2217.5+5943) 
using z = 0.1. 



For the sources where we successfully modeled the syn- 
chrotron spectrum ^ci > ^re- All of these sources show steepen- 
ing of the radio spectrum at higher frequencies which is expected 
in case of spectral aging. For 24P73 the spectral index modeling 
is consistent with the source being a radio phoenix with a total 
source age of 3 x 10^ yrs. We find that the derived synchrotron 
age does not critically depend on the adopted redshift of 0.1 for 
this source (see Table [21). 



6. Discussion 

Based on their toroidal or filamentary morphologies and curved 
radio spectra two of our relics could be classified as ra- 
dio phoenices: VLSS J1515.1-h0424 and VLSS J2217.5-h5943. 
Contrary, VLSS J1133.7-h2324 and the twisted structure in 
IMACS J0717.5-h3745 are probably the result of DSA from struc- 
ture formation shocks. In the case of VLSS J0717.5-h3745, the 
distorted morphology and high temperature of the ICIVL are clear 



evidence for a cluster undergoing a merger. For the other relics 
more observations are needed to determine their origin. 

IMost of the sources in our sample are associated with galax- 
ies in clusters and show irregular morphologies. This can be 
caused by the interaction of the radio plasma with the ICM. The 
steep spectrum of the sources is caused by spectral aging, e.g., in 
some cases the central AGN activity may have stopped, produc- 
ing a so called "dying" radio source. Confinement of the radio 
plasma by the ICM could also have contributed to the steep radio 
spectra, as in this case the radio plasma has "the time" to display 
the eff'ects of spectral aging. 

We have calculated the physical size (largest extent), pro- 
jected distance from the cluster center (7?projected), and 1.4 GHz 
radio power (Pi a) for the radio relics in our sample. We have 
complemented this with values from the literature for radio 
relics with measured spectral indices: A13, A85, A133, and 
A4038 (Sleeetal. 2001), A1240 and A2345 (iBonafede et al " 
2009b), A3667 (Rottgering et al] Il997l) . A548b (Feretti et al 
.2006 ). A2256 (Clarke & Ensslinil2006h . A521 dGiacintucci et al 



12 R. J. van Weeren et aL: Observations of diffuse steep- spectrum sources 

Table 3. Radio relic properties 



Source Name 


1 A/\/\\AlJn, 

^14UUMHz 
"^74MHz 


spectral curvature 

^610MHz ^ 1400MHz 
"^74MHz "^610MHz 


PlA 

10^4 wHz-1 


size^ 
kpc 


-^projected 

kpc 


VLSSJl 133.7+2324 


-1.69 


0.96 


24.1 


570 




VLSSJ1431. 8+1331^ 








125 


175 


24P73 (VLSS J2217.5+5943) 


-2.20 


2.0 








MACS J0717.5+3745 


-1.09^ 


-0.02 


138^ 


700 


320 


Abell 2048 (VLSS J15 15. 1+0424) 


-1.50 


1.60 


0.566 


310 


310 


Abell 2256 


-1.2 




3.95 


1100 


300 


Abell 1240-N 


-1.2 




0.427 


650 


700 


Abell 1240-S 


-1.3 




0.730 


1250 


1100 


Abell 2345-E 


-1.3 




2.62 


1500 


890 


Abell 2345-W 


-1.5 




2.83 


1150 


1000 


Abell 13 


-1.79 


0.65 


0.853 


260 


170 


Abell 85 


-2.30 


1.36 


0.322 


150 


430 


Abell 133 


-1.70 


1.39 


1.07 


55 


40 


Abell 4038 


-1.88 


1.03 


0.0983 


55 


35 


Abell 3667 


-1.1 




17.4 


1920 


1950 


Abell 548b A 


~ -2 




0.260 


310 


500 


Abell 548b B 


< -2.0 




0.250 


370 


430 


Abell 521 


-1.48 




2.90 


1000 


930 


Coma cluster^ 


-1.18 




0.284 


850 


1940 


Abell 2163 


-1.02 




2.23 


450 


1550 


Abell 2744 


-1.1 




6.20 


1620 


1560 


Abell S753 


-2.0 




0.205 


350 


410 


Abell 115 


-1.1 




16.7 


1960 


1510 


Abell 610 


-1.4 




0.444 


330 


310 


RXCJ1314.4-2515-E 


-1.41 




2.08 


920 


685 


RXCJ1314.4-2515-W 


-1.40 




4.81 


920 


685 



^ largest linear scale 

^ excluding the flux of the halo and central head-tail source (see also Ivan Weeren et al.ll2009bl : lBonafede et al.ll2009ah 
' 1253+275 

^ blended with a nearby AGN (except in the GMRT 610 MHz image), therefore no spectral index or spectral curvature could be calculated 



2008 ). 1253+275 ('Giovannini et al."l99ll). A2163 (iFeretti et al 



2004), A2744 (.Orru et al.. .20071 AS753 (ISubrahmanva n et al 



2003 ). A115 (lGovonietal .l l2q0lh A610 dG iovannini & Feretti 
'2000) and RXC J1314.4-2515 (IVenturi et all 2007). The spectral 
indices for our newly found relics were calculated between 74 
and 1400 MHz. The spectral indices for the relics taken from the 
literature were usually measured between 325 and 1400 MHz, 
but for some relics the frequency range is somewhat differ- 
ent (for more details the reader is referred the references given 
above). 

The spectral index of the radio relics, versus the physical 
size is shown in Fig. [20l The projected distance from the clus- 
ter center is color coded. We find that on average the smaller 
relics have steeper spectra. There are two diff'erent explana- 
tions for this correlation. The found correlation between phys- 
ical size and spectral index may (partly) be the result of the pro- 
posed radio phoenices occupying the lower left region in Fig.[20l 
Ignoring the proposed radio phoenices the trend remains, al- 
though the scatter is large so this result is less significant. Such 
a trend though, is in line with predictions from shock statis- 
tics derived from cosmological simulations (lHoeftetal.1 120081: 
ISkillmaneF al. 2008). They find that larger shock waves occur 
mainly in lower-density regions and have larger Mach numbers, 
and consequently shallower spectra. Conversely, smaller shock 
waves are more likely to be found in cluster centers and have 
lower Mach numbers and steeper spectra. We note that more 
spectral index measurements of high quality are needed to con- 
firm the correlation between physical size and spectral index 
(leaving out the proposed radio phoenices). 



The size of the relics versus /^projected is shown in Fig. [2T] 
The projected distance from the cluster center for the shallow 
spectra relics is indeed on average larger than for the steep spec- 
trum relics, indicating that they are mostly located in the pe- 
riphery of clusters where the Mach numbers of the shock are 
higher. We contribute this to the fact that in the periphery the 
shock surfaces are larger and the density/pressure of the ICM is 
lower. Therefore large radio relics are mainly located in the clus- 
ter periphery. Because we are using the projected distance from 
the cluster center, and not the (unknown) de-projected distance, 
some additional scatter is introduced in the plot. Furthermore, 
the spectral indices for the various radio relics are calculated 
using data from radio telescopes with diff'erent array configura- 
tions, sensitivities, and/or frequencies ranges, adding to the scat- 
ter. 

Giant Mpc-scale radio relics are probably caused by DSA 
in an outwards traveling shock front. It is unlikely that they are 
the result of compression and reignition of fossil radio plasma as 
the time to compress such a large radio "ghost" would remove 
most of the electro ns responsible for the rad io emission by radia- 
tive energy losses (IClarke & Ensslinl l2006V In the case of giant 
peripheral relics shock- acceleration is ongoing resulting in rel- 
atively flat spectral indices of about -1, i.e., a balance between 
electron cooling in the post-shock regions and continuous accel- 
eration at the shock front. Behind the shock front, the spectral 
index is indeed observed to steepen for some giant relics (e.g., 
Rottgering et al. 1997). After a few times 10^ yrs the electrons 
behind the shock front have lost most of their energy and cause 
little synchrotron emission. It takes of the order of 1 Gyr for a 
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Fig. 20. Spectral index of radio relics versus their size. Squares 
are the proposed radio phoenices, i.e., the sources from 
ISlee et al.1 (|2001), 24P73 and the relic in Abell 2048. Diamonds 
represent the radio relics tracing merger shocks where particles 
are being accelerated by the DSA mechanism. The color cod- 
ing is according to the projected distance from the cluster center. 
For the relics represented by black symbols we could not obtain 
a reliable projected distance to the cluster center. 



shock wave to travel from the center of the cluster to about the 
virial radius. 

As has been mentioned, an alternative explanation for the 
found size-spectral index correlation could be a possible dif- 
ferent origin of the smaller radio relics. With the smaller radio 
relics originating from the adiabatic compression of fossil radio 
plasma. In fossil radio plasma, the high-frequency synchrotron 
emitting electrons have lost most of their energy. Due to com- 
pression and the increase in magnetic field strength the radio 
plasma becomes visible again. Therefore these sources are char- 
acterized by (very) steep and curved spectra. Proposed example s 
of such sources are 24P73 and those found by lSlee et~al ] (l200Th . 
If spectral aging occurs, the steeper radio spectra should be more 
curved. 

It would be interesting to determine the injection spectral in- 
dices ( Qfinj) for the radio relics. This requires reliable flux mea- 
surements over a wide range of frequencies, especially below the 
break frequency Vb. With enough flux measurements it should 
be possible to separate the efl'ect of spectral aging and a steep 
Qfinj . In particular, if the correlation is explained by the difl'er- 
ent Mach numbers of shocks the injection spectral indices for 
smaller relics located close to the cluster center should be steeper 
than that of larger relics. If this is indeed the case this is evidence 
for DSA. For radio phoenices the injection spectral index should 
be around -0.5 to -0.75 because the fossil radio plasma origi- 
nated from AGN. Their steep spectral indices should then solely 
be the result of spectral aging and not a steep injection spectral 
index. Currently, with a limited number of flux measurements, 
mainly at frequencies where electrons have already lost some of 
their energy, there is a degeneracy between the injection spectral 
index and the amount of spectral aging. New low-frequency ra- 
dio facilities, such as LOFAR, will be needed to determine 
for radio relics. 

The spectral index versus the 1 .4 GHz radio power for the 
relics is shown in Fig.[22l From the figure we can see that there 
is a lack of high power sources with steep radio spectra. This 
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Fig. 21. Projected distance from the cluster center of radio relics 
versus their size. Symbols are defined in Fig.[2Ql 
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Fig. 22. Radio power at 1 .4 GHz of radio relics versus their spec- 
tral index. Symbols are defined in Fig.[20l 



is not surprising since relics with less steep spectra are larger 
(Fig.[20| and therefore should have a higher radio power. 

We have plotted the amount of spectral curvature for our 
sources versus their spectral index, between 74 and 1400 MHz, 
in Fig. [23l The curvature is defined as 0^74-610 - <^6io-i400- A 
higher positive value for the curvature implies that the spectra 
are steeper at higher frequencies. We have only added the four 
sources from Slee et al.l ([2001) for which flux measurements at 
both 1.4 GHz and 74 MHz (VLSS) are available. We calculated 
the 610 MHz fluxes by interpolating the fluxes (linearly in log- 
log space) between 408 and 834 MHz. For other radio relics, 
flux measurements are less reliable or not available at 74 and/or 
610 MHz. As shown in Fig. [23] there is an indication that the 
curvature is higher for the sources with a steeper spectral index, 
although the number of sources in the plot is rather low. If the 
spectral index of relics steepens by spectral aging such a trend 
is expected, but more measurements of spectral curvature will 
be needed to show that this correlation with the spectral index 
really exists. 
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Fig. 23. Curvature of the radio spectra versus their spectral in- 
dex between 74 and 1400 MHz. The curvature is defined as 
^3^74-610 - 0^610-1400- Symbols are defined in Fig.l2Ql 



7. Conclusions 

We carried out 610 MHz radio continuum observations with the 
GMRT of 26 diff'use steep spectrum sources. All of these sources 
were resolved out by 1.4 GHz VLA B-array snapshot observa- 
tions or in the 1.4 GHz FIRST survey. Of these 26 sources 25 
were detected with the GMRT. The radio observations show a 
wide range of sources: radio relics, a giant radio halo, a mini- 
halo, FR-I, head-tail, and USS sources. Here we shortly summa- 
rize the properties of the relics and halos in our sample. 

-The radio relic 24P73 (VLSS J2217.5+5943) is located 
close to the galactic plane. This source was previously known 
to have a steep spectrum but remained unclassified due to the 
limited sensitivity and resolution of previous observations. Due 
to the crowed Milky Way field and high extinction no cluster 
is found in POSS-II images. Deep NIR imaging will be needed 
to reveal the presence of a cluster. The sources has an extremely 
steep curved spectrum, likely the result of adiabatic compression 
of fossil radio plasma and spectral aging. 

-VLSS Jl 133.7+2324 is a filamentary relic located next to a 
nearby spiral galaxy. The radio source breaks up into two par- 
allel filamentary structures. This relic source is probably associ- 
ated with a structure of galaxies located at a distance of z ~ 0.6. 
The nearby spiral galaxy itself is also detected in the radio im- 
ages, the flux density consistent with that predicated by the FIR- 
radio correlation (using IRAS fluxes for this galaxy). The galaxy 
is partly blocking our view of the distant structure of galaxies, 
making follow-up observations more difl&cult. 

-VLSS Jl 43 1.8+1331 consists of two components. The east- 
ern part is probably associated with a cD galaxy in the center of 
the cluster MaxBCG J217.95869+13.53470 at z = 0.16. The 
radio morphology suggests we are seeing signs of interaction 
between the radio plasma and the ICM. The other part is proba- 
bly a relic source. These sources makes an interesting target for 
follow-up X-ray observations. 

-VLSS J0004.9-3457 is associated with an elliptical galaxy 
in a small cluster at z ~ 0.3. The difl'use radio emission sur- 
rounds the galaxy in the form of a mini-halo. The source has an 
arc-like extension to the east. 



- VLSS J07 17.5+3745 is located in the cluster 
MACS J07 17.5+3745, a massive merging system at a red- 
shift of 0.5548. The radio emission is complex, consisting of a 
large elongated radio relic and a giant 1.2 Mpc radio halo. The 
halo has a radio power (Pia) of 5 x 10^^ W Hz"\ the highest 
one known to date. The relic might trace a giant shock wave 
related to the cluster merger. This is consistent with the very hot 
ICM and relatively flat spectral index. 

-VLSS J15 15. 1+0424 is a radio relic consisting of three 
filamentary structures in the periphery of the cluster Abell 2048. 
Adiabatic compression of fossil radio plasma probably resulted 
in the complex morphology. 

We find that larger relics are mostly located in the cluster 
periphery, while smaller relics are found closer to the cluster 
center. We also discovered a correlation between the spectral in- 
dex and the physical size of the relics. A likely explanation for 
this correlation is that the larger shock waves occur mainly in 
lower-density regions and have larger Mach numbers. As a con- 
sequence they have shallower spectra. However, the correlation 
can also (partly) be explained by invoking diff'erent origins for 
relics. Relics formed by diff'usive shock acceleration are direct 
tracers of large shock fronts in clusters, while the compression of 
fossil radio plasma is thought to produce relatively small relics 
with steep curved spectra. 

If the correlation is explained by the fact that the larger shock 
waves occur mainly in lower-density regions and have larger 
Mach numbers then the injection spectral indices should flatten 
away from the cluster center for the various relics. To measure 
this requires a sufficient number of flux measurements, espe- 
cially below the break frequency in the spectrum. With upcom- 
ing new radio facilities operating at lower frequencies it should 
therefore be possible to break the degeneracy between the injec- 
tion spectral index and spectral aging. 

Follow-up observations to measure the polarization proper- 
ties are currently underway. We also plan to create spectral in- 
dex maps for our sources which should give more information 
on the life-times of the various synchrotron emitting regions. 
Optical observations have been taken to characterize the large- 
scale environment around the radio sources. For two of our dif- 
fuse sources these observations should confirm the presence of 
a cluster. X-ray observations will be needed to study the dynam- 
ical state of the clusters and reveal the presence of shock fronts 
in the ICM. 
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Appendix A: Other sources in the sample 

A.1. VLSSJ0227.4-1642 

The radio map shows two resolved elongated patches of emis- 
sion. The brighter eastern component has a flux of 25.8 mJy 
and the weaker western component 1 1.2 mJy. Both components 
are resolved, 28.3'' by 10.3'' (eastern component) and 22.3" by 
9.0" (western component), and show structure at the limit of the 
~ 5" beam. The position angles of the longest axes are roughly 
the same, which suggests a link between the two components. 
No optical counterpart is detected at the limit of the POSS-II sur- 
vey. The steep spectral index of about -1.5 could be explained 
if this source is a FR-I type AGN, where the central engine has 
stopped and the synchrotron emission in the two lobes has steep- 
ened by spectral aging. In this case both components should 
roughly have the same spectral index. 



A.2. VLSSJ0250.5-1247 

The radio observations show a double source. The two compo- 
nents are slightly resolved, have a similar flux and seem to be 
connected by a faint bridge. The separation between the two ra- 
dio components is 15". No optical counterpart is detected. 



A.3. VLSS J0646.8-0722 

The radio map shows a source with an angular size of 1.0' by 
0.5'. An optical counterpart (PCG 76039, mag/?=18.4) is de- 
tected in both the POSS-II and 2MASS surveys. This galaxy 
was also listed in the catalog of Galaxi es Behind the Milky 
Way (CGMW 1-0379, ISaito etal.1 11990 ') and in the catalog 
Galaxies in the "zone of avoidance" (ZOAG G218.99-04.37, 
[Weinberger et al.lll995[i) and was in both cases classified as an 
elliptical galaxy. The amount of extinction in this region is 
Afi = 2.2 mag. The radio map shows a one-sided tail originating 
from the optical counterpart. This might indicate that the source 
is located in a cluster and we see the eff'ects of the interaction 
between radio plasma and the surrounding ICM. 



A.4. VLSS J111 7.1 +7003 

This source has a remarkable steep and straight radio spectrum 
without any indication of a flux turnover at lower frequencies. 
The source is resolved into a smooth featureless roughly spheri- 
cal blob ( 26" by 23"). No optical counterpart is detected. Given 
the steep straight radio spectrum it could be a radio halo, how- 
ever the small size of the source argues against such an iden- 
tification. The source could be a mini-radio halo surrounding a 
central cD galaxy of a distant (proto) cluster. 



A.5. VLSS J1 636.5+3326 

The radio map shows two compact components, with the north- 
ern component extending to the west. SDSS images show several 
faint mostly red galaxies around the radio source. Photometric 
redshifts (SDSS DR7) indicate the presence of a cluster at z = 
0.65 ± 0.19. Several galaxies can be identified as possible optical 
counterparts. 



A.6. VLSS J1 71 0.5+6844 



This "head-tail" radio source has a length of 0.73' and a width 
of 0.27'. An optical counterpart is detected, a galaxy with a K- 
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magnitude of 14.75. An optical /-band image from the 2.5m 
Isaac Newton Telescope (INT) shows that this source is located 
within a cluster at a redshift of z ~ 0.3. An east- west elon- 
gated X-ray source IRXS J17 1034.0+684403 follows roughly 
the galaxy distribution confirming the presence of a cluster. 

M. VLSSJ1 930.4+1 048 

The radio map shows an elongated source with a LAS of lOT'. 
The radio morphology is consistent with a FR-I radio source. No 
optical counterpart is detected at the position of the core, but this 
is not surprising given the amount of extinction (Ab = 3.759) at 
this galactic latitude of -3.6 degrees and the high star density 
in the POSS-II images. A 69.6 mJy (probably unrelated) radio 
source is located 70''to the north. 



A.8. VLSSJ2043.9-1118 

An optical counterpart for this source is detected in POSS-II im- 
ages. The faint optical source is located on the peak of the radio 
emission. The radio emission surrounds the galaxy and there is 
a hint of a faint extension to the east. 



A.9. VLSS J2044.7+0447 

The source is elongated in the east-west direction and has a 
LAS of lOr'. A R-magnitude 19.56 galaxy corresponds to a 
peak in the radio emission, consistent with this being a FR-I 
radio source. The high rms noise is the radio map is caused 
by dynamic range limitations from 4C 04.71, a 1.24 Jy source 
(610 MHz), located 3.3' to the south. 

A.10. VLSS 2122.9+0012 

The radio source is compact, 11.6'' by 10.7". This is probably 
a USS source associated with a HzRG. The lack of an optical 
counterpart in SDSS images in consistent with this explanation, 
and implies the source is located at z > 1.4. 



A.I 3. VLSSJ2345.2+2157 

The radio source shows an S- shape symmetry. An optical coun- 
terpart is detected in 2MASS and POSS-II images being the cen- 
tral cD galaxy of a cluster. This is confirmed by the presence of 
the X-ray source IRXS J234518.4-h215753. 

A.I 4. VLSS J051 1.6+0254 

The radio map shows a FR-I type radio source. Two possible 
optical counterparts are detected which are the central elliptical 
galaxies in a cluster at a redshift of about 0.2. 

A.I 5. VLSS J051 6.2+01 03 

VLSS J05 16.2-^0103 is a relatively compact source, without a 
double lobe structure. The source does not have an optical coun- 
terpart in POSS-II images. This implies the source is probably 
located at z > 0.7. 

A.I 6. VLSS J2209.5+1 546 

The radio map shows an elongated source. No optical counter- 
part is detected at the limit of the POSS-II images. The source is 
probably a distant FR-I source. Optical imaging will be needed 
to confirm this classification. 

A.I 7. VLSSJ2241.3-1626 

This source could be a disturbed FR-I source, no optical counter- 
part is detected. Optical imaging will be needed to confirm this 
classification. 

A.I 8. VLSS J2357. 0+0441 

The source is found to be a blend of three compact sources, all 
without optical counterparts. By convolving the GMRT image 
to the NVSS resolution (45") it seems that the two outer compo- 
nents are responsible for the steep spectrum. These are probably 
distant FR-II sources. 



A.11. VLSS J221 3.2+3411 

This source vs^ a s als o included in the USS sample of 
'De Breuck et al." ('2000), WN J2213+3411. A K, band image 
from De Breuck et al. ( 2002b) showed a possible counterpart. 
However, because no high-resolution (< 10") VLA images were 
available the identification was uncertain on the basis of position 
alone. Our high-resolution 5" image shows a distorted double- 
lobe source. The central radio core is not detected. The NIR 
counterpart is located in between the two lobes, confirming that 
this is indeed the host galaxy. The K-band magnitude of 18.3 im- 
plies a redshift of 1 .6 ± 0.5 for this source using the K-z relation. 



A.I 2. VLSS J2341. 1 + 1 231 

The radio map shows a FR-I type radio source. The northwest- 
ern and southeastern lobes are Hsted as NVSS J234104-h 123203 
and NVSS J234107+123126, respectively. An optical counter- 
part is detected on POSS-II images at roughly the expected po- 
sition based on the location of the lobes. This source may be a 
so called "dying" radio source. 
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Fig. A.l. GMRT 610 MHz radio map. Contour levels are drawn as in Fig.[T] 
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Fig. A.2. GMRT 610 MHz radio map. Contour levels are drawn as in Fig.[T] 
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Fig.A.3. Optical images for radio sources with counterparts. Images are from SDSS or POSS-II (red), except for 
VLSS J17 10.5+6844 for which we took seven /-band images from the INT Wide Field Camera with a total exposure time of 
2520 sec. GMRT 610 MHz contour levels are drawn as in Fig. [21 



